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Abstract 

Trial wavefunctions, constructed explicitly from the unique 2-dimensional 
Mott insulating state with antiferromagnetic order, are proposed to describe 
the low-energy states of a Mott insulator slightly doped with holes or elec- 
trons. With the state behaving like charged quasi-particles with well-defined 
momenta, a rigid band is observed. These states have much less pairing cor- 
relations than previously studied ones. Small Fermi patches obtained are 
consistent with recent experiments on high T c cuprates doped lightly with 
holes or electrons. States showing the incoherent and spin-bag behaviors are 
also discussed. Using these wavefunctions, a number of results obtained by 

exact calculations are reproduced. 
PACS: 74.72.Jt, 75.50.Ee, 79.60.-i 
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With the continuous improvement of experimental technique and sample making, it 
has recently become possible to study in detail the phenomena in the lightly doped high 
T c cuprates. Indeed, many intriguing behaviors concerning the physics of doping the 
two-dimensional (2D) Mott insulator in this very undedoped regime are observed. Re- 
cent angle-resolved photoemission spectroscopy (ARPES) results unearth contrasting be- 
haviors between lightly hole-doped C a<i- z N a z CuO<iCl<i (iVa-CCOC) [p]] and electron-doped 
Nd 2 - x Ce x Cu0 2 (NCCO) high T c cuprates @. Although ARPES on the undoped (i.e. 
z=x=0) insulating state shows an identical energy dispersion of a single hole created below 
the charge gap, results at a little higher dopings are demonstrated to be different: while 
a small hole patch is observed to be centered clearly at momentum (7r/2,7r/2) in the Na- 
CCOC at even z=0.1, small electron patches centered at (vr,0) and (0,7r) are observed for 
£=0.04 in NCCO. Moreover, the in-plane transport in lightly doped systems shows high 
mobility of the charge carriers j3]. This and the recent high- resolution scanning tunneling 
microscopy/spectroscopy results indicate the signature of the quasi-particle behavior Q. 

The single hole/electron behavior and its dispersion have been studied using various 
approaches on the t-t'-t"-J model P,|6|,|7|,|8|J^JT(]|] . However, these studies on doping holes and 
electrons into the system only emphasize the asymmetry resulting from the different signs 
of t' and t" for the corresponding Hamiltonians. It is unclear whether the same physics is 
working for these two systems with different Hamiltonians. Should one try to construct a 
different theory when electron-doped cuprates are considered? Furthermore, do the models 
predict small Fermi surfaces and quasi-particles? 

There are other issues concerning the model itself which have been bothersome. The stud- 
ies applying the projected d-wave superconducting, or the resonating- valence-bond (RVB), 
state with the antiferromagnetic long range order (AFLRO) also suggest that away from 
half-filling, the superconductivity (SC) revives and the ground state shows both SC and 



AFLRO [T^|TJ]. But so far most experiments do not support the coexistence. In addition, 
there has been not enough understanding on various properties found in the exact results of 
the t-J type models on finite clusters |]5l,|6|J^j9|,p^, |rT| , p^1 in terms of the many-particle wave 



functions (WF's). 

In this paper, we discuss our recent proposal of a theory based on the variational approach 
to understand these issues [p~5] , p~6|] . Specific trial wavef unctions (TWF's) are constructed to 



describe both the low-energy states of the associated t-t'-t"-J models with lightly doped holes 
and electrons. These WF's are generalizations of the single-hole WF first written down by 
Lee and Shih [0. In contrast to other TWF's [fTTH , ours are constructed solely from the 
optimized one at half- filling and include no hopping amplitudes t' and t". Yet, surprisingly, 
a number of properties including dispersion relations, momemtum distributions, spectral 
weight etc. are obtained correctly for both hole-doped and electron-doped systems. In the 
following, the Lee-Shih WF for one doped hole is used and generalized for several holes as 
well as electron-doped systems. All numerical results reported below are from variational 
Monte Carlo (VMC) calculations for an 8 x 8 lattice with periodic boundary conditions. 
At half-filling, the system is described by the Heisenberg Hamiltonian 

Uj = J^(Si • Sj - -niUj), (1) 

where denotes nearest-neighbor (n.n.) sites. Each site is occupied by only one sin- 

gle electron. As holes or electrons are doped into the system, we consider the Hamilto- 
nian TC=TC t -t>-t"+7~Lj including n.n. and longer range hoppings. Here 7~lt-t'-t" = —tY,(i,j)a 
cl a Cj a - t'Y,{i,i)o c\ a ci a - t"J2(i, m )cr c| <T c mo -+H.c. with and (i,m) representing the sec- 
ond n.n. and third n.n. site pairs. Note that Cj CT in 7i t _ t i_ t n creates different kind of holes 
from single-electron-occupied sites at half-filling: empty holes (Oe-hole) for hole doping and 
two-electron-occupied holes (2e-hole) for electron doping flSfl . Operator c ia is actually equal 
to c ia (l — rij -o-) or Ci_ a n ia for hole or electron doped case, respectively. Therefore, despite 
the constraints, states in the two cases are in one-to-one correspondence after a local trans- 
formation Ci a — > c|_ CT is made. However, because of the Fermi statistics, the exchange of a 
single spin with a 2 e- hole has an extra minus sign as compared to the Oe-hole. Hence, the 
only difference between the hole-doped and electron-doped t-t'-t"-J model is t'/t — > —t'/t 
and t" jt — ► —t"jt after we change the Cj CT on B sublattice sites to — Cj CT 0. With all these, 



we then treat the hole and electron doped cases in the same manner with the Lee-Shih WF 
originally proposed only for a single hole. The VMC results presented in this paper are for 
J/t = 0.3, if ji = — (+)0.3 and if' ji = +(— )0.2 in the hole (electron) doped case following 
the values usually used B. 



We shall apply the standard VMC method [T7] that enforces the local constraint exactly. 
Following the work by Lee and Shih , the TWF for one doped hole with momentum q and 
S z =l/2 is constructed to have {N s /2) — 1 singlet pairs of electrons and a single unpaired 
electron with momentum q and S z =l/2, 

k 

The prime on the summation symbol indicates that the momentum q is excluded from the 
sum if q is within the sublattice Brillouin zone (SBZ), otherwise, q — Q is excluded. N s 
here is the total number of sites and Q=(vr, ir). 

is explicitly constructed from the optimized half-filled WF |^ ) = 
-PdEk(^k a kt a -k| + -^k^kT^-k|)] Ars//2 |0) an d h does not contain any information about hop- 
pings, t', t" and neither explicitly t, of the doped hole or electron. However, the effect of t is 
included in the RVB uniform bond X=(J2a c la c jcr) which describes the large quantum fluc- 
tuation and spin singlet formation. There is also no need to introduce t' and t" in the TWF 
as they are compatible with AFLRO. Here the coefficients and are functions of £k 
andA k . ±£ k = ±(e k + ( Jm s f) 2 are energy dispersions for the two spin density wave (SDW) 
bands with e^=— jJx(cosk x + cosk y ) and the staggered magnetization m s =(S^)=— (Sg), 
where the lattice is divided into A and B sublattices. a^ a and represent the operators 
of the lower and upper SDW bands, respectively, and are related to the original electron 
operators and Ck+Qo- with Q=(7r,7r) set for the commensurate SDW state. A k = |JAd k 
with =cosk x — cosk y here is for the (i-wave RVB (d-RVB) order parameter. The projection 
operator Pd enforces the constraint of no doubly occupied (or vacant) sites for cases with 
finite hole (or electron) doping . At half-filling, iV s equals the total number of electrons. 
Notice that the sum in |^ ) is taken over SBZ. There are two variational parameters: A/x 



and m s /x in these WF's. 

The energy dispersion obtained from for one doped hole has been shown by Lee and 
Shih |7] to agree very well with that of exact calculations, self-consistent Born approxima- 
tions (SCBA), and Green function Monte Carlo methods etc.. As for the case of having an 
extra up-spin electron with momentum q doped into the half-filled state, the energy disper- 
sion can be calculated with this same WF The only difference is signs of t'/t and t"/t 
are changed in the Hamiltonian. 

The variational energies for one doped electron are shown as black dots in Fig.|l| This 
result agrees well with that of Xiang and Wheatley |J obtained by SCBA. The optimal 
variational parameters are (A/x, m s /x) = (0- 25, 0.125). The ground state is at momentum 
q = (vr,0). The VMC results can be fitted simply by Ei k =E k — 2t e ff(cosk x + cosk y ) — 
4t' e jj-cosk x cosk y — 2t" fj (cos (2k x ) + cos(2k y )) with parameters described in the caption of 
Fig.|l]. The dispersion thus seems to be simply the combination of the mean-field band at 
half-filling and the coherent hoppings [7J]. 

To examine further the physical properties of we calculated the momentum distri- 
bution function (n^(k)) for the ground state of a single hole with momentum q = (tt/2, 7r/2) 
and S z =l/2. Results are shown in Fig.|2|(a) and (b). Note that the dips or pockets at 
(7r/2, 7r/2) and anti-dips at (— tt/2, — vr/2) found by Leung for the exact results of 32 
sites are also clearly seen here. It is quite amazing that including no if and t", not only 
produces the correct energy dispersions for a single doped hole or electron it also provides a 
correct picture about the momentum distribution. 

The momentum distribution functions (n^(k)) for electron doped systems could be also 
calculated from if we perform the transformation, Cj CT — > c\_ a , and t is chosen to be 
positive. In fact, it is easy to show that (n® (k)) = 1 — (r^ CT (Q — k)). The results for the 
ground state of a single doped electron with momentum k = (vr, 0) and spin S z =l/2 are 
shown in Fig.|2|(c) and (d). There are peaks at k = (n, 0) and an anti-peak at k = (0, ir). 

Now we shall generalize the Lee-Shih WF j^i) to the case of two holes. The simplest 
possible way is just to take out the unpaired spin from if we are interested in the state 
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with zero total momentum and S z =0, which turns out to be the lowest energy state. The 
TWF for two holes with momenta q and — q is 

k 

Note that the momentum q is not included in the summation. It is most surprising to find 
that although has zero total momentum irrespective of q, its energy varies with the 
missing momentum or the hole momentum q. The dispersion turns out to be very similar 
to that of a single electron as shown in Figjj] ||15|| . The state with momentum q=(7r, 0) has 
the lowest energy for two electrons. The values of the two parameters A/x and m s /x are 
the same for |^ 2 ) an d 

Using WF 1^2); the energy dispersion for two holes doped into the half-filled state again 
has an almost identical form as that of a single hole and the minimum is at q = (vr/2, 7r/2) 
0. The lowest energy obtained is — 26.438(3)t which is much lower than the variational 



energy, — 25.72(l)i, using the TWF applied by Himeda and Ogata jTJj]. Even if we include t' 



and t", the variational energy — 25.763(7)t, is still much higher than ours [H5|. We also found 



that the hopping amplitudes for n.n., second n.n. and third n.n. of two holes are almost 
twice that of one hole. The momentum distribution function for this state (not shown) has 
dips at (7r/2, 7r/2) and (— 7r/2, —tt/2). This is in good agreement with the exact result [Hj 
for the t-t'-t"-J model with 2 holes in 32 sites. 

It is then straightforward to write down the same type of TWF for three holes with 
momenta q, q' and — q': 

1*3) = p d 4 T [E'(^4Ai + 5 k &k T ^)] (iVs/2) - 2 |o), 

k 

where q' and q are excluded from the summation. Just like the case with two holes or two 
electrons energy dispersions are proportional to the sum of the three single hole energies 
at momenta q, q' and — q'. In FigJ^, the dispersions are plotted as functions of q. For 
the electron doped {i.e. with three 2e-holes) case in Fig.[5](a), the lowest energy state is at 
q = (37r/4, 0) (shown as an open circle) within the SBZ after the first two electrons occupy 
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q' = (tt, 0). With three doped holes (FigfJ(b)), after the two electrons at (tt/2, tt/2) are 
removed the ground state is now at (tt/2, — tt/2) (the open circle). As shown clearly in 
Fig.|3], the dispersions follow nicely the single hole one. 

Values of the staggered magnetization m=N s ~ 1 J2i(— l)*^* are computed and compared 
for several Oe-hole and 2e-hole concentrations. With the same variational parameters, it 
is found that the preference of (tt/2, tt/2) for Oe-holes causes clearly larger disturbance of 
the AF order than for the electron doped case where 2e-holes with momentum (it, 0) shows 
much less influence on the AF order ||15|| . This is consistent with previous work ||. It is 
also consistent with experimental results that AF phase is more stable for electron doping 
than hole doping | p0[| . 

So far, based on the t-t'-t"-J model we have proposed a TWF to describe the low energy 
states of slightly doped antiferromagnetic Mott insulators. Exactly the same TWF's are 
proposed to account for the behavior of both hole doping and electron doping, after we 
employed the particle-hole transformation. Different energy dispersions for these two cases 
are due to the different signs of t' /t and t"/t which is a direct consequence of the constraint 
that electron doped system has 2e-holes while hole doped system only has Oe-holes. Rigid 
band and quasi-particle behavior are demonstrated for both cases. The theory provides an 
explanation of recent ARPES results. In lightly hole-doped cuprates, small Fermi pocket is 
around (tt/2, tt/2). In electron-doped cuprates, the patch is around (tt, 0) as shown in the 
inset of Fig.[l| with doping at about 3%. 

Another important property of our WF's is that holes are essentially independent of 
each other as they obey the same energy dispersions (with very little renormalization of 
parameters). Exactly because this quasi-particle like property is unchanged after doping, 
our WF has AFLRO but very little superconducting pairing correlations. The presence of 
superconducting state certainly will change the excitation spectra. In particular, the ci-wave 
SC, which coexists with AFLRO in some of the previous variational studies, should have low 
energy excitations along the nodes. This is certainly not seen in our TWF's. In addition, 
the holes are not attractive to each other in our WF's. In Fig.f| the hole-hole correlation 



function for our TWF and the WF used by authors in |14| are compared for 4 holes in 
an 8x8 lattice. The lack of attraction between holes is consistent with Leung's low energy 
states obtained exactly for two holes in 32 sites |I0[ . Long range (i-wave pairing correlation 
defined ]T7|] for our TWF and that in Himeda-Ogata one are, on average, about 0.002, and 
0.018, respectively. Thus the c?-RVB pairing for spins assumed by our WF's are not in any 
way implying the pairing of charges. 

Although our TWF's or the Lee-Shih WF's has reproduced many numerical results ob- 
tained by exact diagonalization, SCBA, etc. for one or two holes or electrons, there exists, 
however, some inconsistancy in the detail of the comparison. Namely, at some momenta 
k's the spectral weights Z k = |(\l f k|ckcr|^ f o)| 2 /(^ol c k - c ko-|^ / o) computed for the one doped 
hole case using |^ k ) = are much larger than the exact results J8|. This may indicate 
that the quasi-particle states do not exist everywhere in the Brillouin zone. In fact, the 
ARPES on the undoped (z = 0) Na-CCOC also show well-defined peaks only locally in the 
high-symmetry directions Jl[. 

We have constructed TWF's with which small spectral weights are realized at particular 
k points [T(J. For the single hole doped case, state described by WF 



l*i> = p* <A E '(^4/- k i + bA^i)] {Ns/2) - 1 \°) 

[Mqft] 

with the hole momentum q/j = (tt/2, 7r/2) has also one singlet bond less than the half-filled 
case, but the unpaired spin momentum q s is chosen here not to be the same as q^. Note 
that only q/j is excluded within the sum for 1^). With the same parameters A/% and m s /xi 
it is found that can have lower variational energy than that of the quasi-particle state 
at some q s 's. Apparently, many states could be constructed with same q s but different 
q^. These states constitutes the continuum at higher energy. 

Substituting WF and the quasi-particle into |^k) and computing Z k suffices 
to reproduce the variation of the spectral weights obtained in the exact results ||. States 
described by 1^^) have smaller Z^s and are thus with a strong incoherent character. We 
also found that more quasi-particle states have lower energy below the continuum in the 
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pure t-J model with t' = t" = JTTJ than in the t-t'-t"-J model. This is consistent with 
what has been known in the exact results ||. 

To understand the incoherent states more, the hole-spin correlations for with dif- 
ferent q s 's are examined [0. The spin configurations around the hole are quite different 
from that of the quasi-particle states. Spin moments around the hole in the quasi-particle 
state can have values larger than the average at that momentum, thus the unpaired 
spin seems to be bound to the hole. However, in the incoherent state the staggered magne- 
tization is suppressed around the hole, this is like a spin-bag state |2lj. The unpaired spin is 
no longer associated with the hole and this is now a spin-charge separated state. The amaz- 
ing difference between these two wave functions could be understood quite easily. When 
we apply the spin flipping operator, S^(\s.'; q' = q[= q h \) = £ q ' c q , +k , m) c q > i(T) | q / =q with 
q + k = q s , to the quasi-particle state, the unpaired spin q is excited to a new momentum 
q s in the state \^f'x). Thus the unpaired spin is no longer bound to the hole. Furthermore, 
the spin-spin correlations across the hole with q s = (tt, 0) show the anti-phase domain, i.e. 
AF correlation at the same sublattice. These features reproduce what have been obtained 
in the exact calculations |3|. 

In summary, we have discussed new TWF's describing the low-energy states of the t-t'- 
t"-J model at lightly doping. These states reproduce various exact numerical results and, 
also, show consistent behaviors with what have been observed in the experiments. 
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FIGURES 

FIG. 1. Energy dispersion of one electron in the t-t'-t"-J model on an 8 x 8 lattice. Black dots 
are VMC results by using The fitted dispersion Eik — Eq are plotted as gray diamonds with 

parameters x = 6.92, A = 2.71, m s = 18.84, E = 7.43, t eff = 0.06, t' eff = -0.15, t" eff = 0.1. 
Inset: patches in one of BZ by filling the fitted dispersion in the main figure up to ~ 3% doping. 

FIG. 2. Momentum distribution functions (rv^(k)) for a single hole, (a) and (b), and electron, 
(c) and (d), in the 8x8 t-t'-t"-J model. A scale is shown in between each set. The darker area 
indicates smaller values of (n^ e \k)}. 

FIG. 3. Energy dispersions of (a) three doped electrons and (b) three doped holes in the 
t-t'-t" -J model. Filled circles are VMC results for the momenta shown along the horizontal axis 
by using A total minus sign has been multiplied to the hole doped case. Open circles, with 

the momenta also displayed, represent the lowest energy states. Results here are obtained with 
parameters shown in the figures. 

FIG. 4. Hole-hole correlation functions with 4 doped holes in the 8x8 lattice. The result 
obtained using our WF, i.e. same form as \^2) with momenta (7r/2,7r/2) and (tt/2,— tt/2) excluded 
in the sum, is compared with that using the Himeda-Ogata WF . Results here are obtained with 
optimized parameters (A/x, m s /x)=(0. 25, 0.125) and, for Himeda-Ogata one, also /i=-0.025. 
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